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ABSTRACT 

Observations of GRB 100724B with the Fermi Gamma-Ray Burst Monitor (GBM) find that the 
spectrum is dominated by the typical Band functional form, which is usually taken to represent a 
non-thermal emission component, but also includes a statistically highly significant thermal spectral 
contribution. The simultaneous observation of the thermal and non-thermal components allows us to 
confidently identify the two emission components. The fact that these seem to vary independently 
favors the idea that the thermal component is of photospheric origin while the dominant non-thermal 
emission occurs at larger radii. Our results imply either a very high efficiency for the non-thermal 
process, or a very small size of the region at the base of the flow, both quite challenging for the 
standard fireball model. These problems are resolved if the jet is initially highly magnetized and has 
a substantial Poynting flux. 

Subject headings: Gamma-ray burst: individual: GRB100724B - Gamma rays: stars - Radiation 
mechanisms: thermal - Radiation mechanisms: non-thermal - Acceleration of par- 
ticles 



1. INTRODUCTION 

The prompt emission detected from Gamma-Ray 
Bursts (GRBs) is believed to originate at large dis- 
tances from the ce ntral engine , from within an ultra- 
relativistic outflow (Piran 2004). This ultra-relativistic 
motion is necessary to avoid strong 77 annihilation, a 

sylvain.guiriec@nasa.gov, sylvain.guiricc@uah.edu 

± University of Alabama in Huntsville, NSSTC, 320 Sparkman 
Drive, Huntsville, AL 35805, USA 

2 Department of Physics, Royal Institute of Technology, Al- 
baNova, SE-106 91 Stockholm, Sweden 

3 The Oskar Klein Centre for Cosmo Particle Physics, Al- 
baNova, SE-106 91 Stockholm, Sweden 

4 Institut dAstrophysique de Paris UMR 7095 Universite 
Pierre et Marie Curie-Paris 06; CNRS 98 bis bd Arago, 75014 
Paris, France 

5 Institut Universitaire de France 

6 Dept. of Astronomy & Astrophysics, Dept. of Physics and 
Center for Particle Astrophysics, Pennsylvania State University, 
University Park, PA 16802, USA 

7 NASA Goddard Space Flight Center, Greenbelt, MD 20771, 
USA 

8 Department of Physics and Department of Astronomy, Uni- 
versity of Maryland, College Park, MD 20742, USA 

9 Istituto Nazionale di Fisica Nucleare, Sezione di Pisa, I- 
56127 Pisa, Italy 

10 Space Science Office, VP62, NASA/Marshall Space Flight 
Center, Huntsville, AL 35812, USA 

11 NASA Postdoctoral Program Fellow, NASA/Marshall 
Space Flight Center, 320 Sparkman Drive, Huntsville, AL 35805 

12 Universities Space Research Association, NSSTC, 320 
Sparkman Drive, Huntsville, AL 35805, USA 

13 Max-Planck-Institut fur extraterrestrische Physik (Giessen- 
bachstrasse 1, 85748 Garching, Germany) 

14 National Space Science and Technology Center, 320 Spark- 
man Drive, Huntsville, AL 35805 

15 Jacobs Technology, Inc., Huntsville, Alabama 

16 Los Alamos National Laboratory, PO Box 1663, Los 
Alamos, NM 87545, USA 

17 University College, Dublin, Bclficld, Stillorgan Road, 
Dublin 4, Ireland 



signature that is not observed (see, e.g., |Pirar] [r999). 
Thermal emission is naturally expected in such a sce- 
nario. Indeed, since the densities at the base of the 
relativistic flow are very large, the medium is optically 
thick to radiation owing to Thomson scattering by en- 
trained electrons. The optical depth decreases during 
the relativistic expansion and the outflow eventually be- 
comes transparent for its own radiation, at the photo- 
spheric radius. Any internal energy that is still car- 
ried out by the flow can be radiated at the photosphere 
and will be observed as a thermal component in the 
prompt spectrum. This expected photospheric emis- 
sion in GRB spectr a was early sugg e sted on such the - 
oret ical grounds bv iGoodmanl (|1986f ). iMeszarosI {2002), 
and lRees and M eszaros (|2005l ). among others. The non- 
thermal component observed in the spectrum has to be 
produced by another mechanism in the optically thin re- 
gion, i.e., well above the photosphere. Due to the ultra- 
relativistic motion, this difference in the radius of the 
emission implies a delay between the observation of the 
two components that is usually small compared to the 
typical duration of a long GRB, and is also small com- 
pared to the typical duration of time intervals used for 
time-dependent spectroscopic analysis. The thermal and 
non-thermal components sho uld then appear superim - 
posed for the observer (e.g. | Meszaros and Reed l2000f ). 
iDaigne and Mochk ovitchl (|2002f ) pointed out that in the 
standard fireball model, the photospheric component can 
easily be dominant in the spectrum if the efficiency /nt 
of the mechanism responsible for the non-thermal emis- 
sion is only moder ate (/nt 40%) ■ 
Ob se^va tionally, [ GhirlandaelaLj (|2003f l. iRvdd (|2004l 
|2005[ ). and lRyde et alJ (|2010f l argued that a photospheric 
component is present in CGRO BATSE data. The lim- 
ited energy range provided by BATSE (20-2000 keV), 
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however, hampered the possibility of unambiguously 
identifying the emission process. Since the launch of 
Fermi in 2008, the combination of the Gamma Ray Burst 
Monitor (GBM) and the Large Area Telescope (LAT) 
provides an unprecedented energy range for GRB spec- 
troscopy, and the identification of the emission processes 
responsible for the gamma-ray prompt emission may be- 
come a reality. GBM alone covers a wider energy range 
than its predecessor BATSE, and the design of its data 
enables finer resolution spectroscopy. This allows better 
constraints on spectral fits, with increasingly complex 
models. 

GRB energy spectra in the keV-MeV energy range 
are usually well re presented by the Band func- 
tion ([Band et al.lll993D . two power-laws, smoothly joined 
and parameterized by E pca k which represents the energy 
at which peak power is radiated (jGehrelsl 119971 ) . The 
value of the low-energy power-law index, a, is higher 
than the value of the high-energy power- law index, j3, 
and the parameter E pca ]< of the Band function for GRBs 
generally app ears to follow predictable trends wit h time 
and flux level (|Ford et al.ll995HGuiriec et al.ll2010f ). It is, 
however, an empirical function rather than a physically- 
motivated model. The meaning of the parameters in 
the context of emission and transport mechanisms taking 
place in GRBs is not well understood, but is generally 
believed to represent the non-thermal emission from ac- 
celerated charged particles. 

In section 1, we describe our observational results con- 
sisting of a GRB prompt-emission spectrum best fit with 
the combination of a thermal component and a standard 
Band function. In section 2, we use these results to con- 
strain the origin of the energy released in the GRB jet. 

2. OBSERVATIONS 

GBM is composed of 12 sodium iodide (Nal) detectors 
covering an energy range from 8 keV to 1 MeV and two 
bismuth germanate ( BGO) detectors sens itive between 
200 keV and 40 MeV (jMeegan et al.ll2009D . The instru- 
ment triggered on 24 July 2010, at Tn=00:42 :05.992 UT 
on the very bright GRB 100724B (|BhatH2010( ). The event 
was also seen at hi gher energies in the Fermi Large Area 
Telescope (LAT) (|Tanaka et al.l [2010h . The most pre- 
cise position for the direction of the burst is the inter- 
section of the Interplaneta ry Network annulus obt ained 
using GBM, Konus-WIND ((Golenetskii et al.ll2010T ). and 
MESSENGER data with the 90% LAT confidence level 
location error box, and is a strip of sky centered on RA = 
118.8° and Dec = 75.8° which is 1.2° long and 0.2° wide 
(K. Hurley and V.Pal'shin, private communication). Fig- 
ure Q] (top two panels) shows the GBM light curve of 
GRB 100724B in two energy bands. Multiple peaks 
of varying intensity are superimposed on a pre-trigger 
plateau, with a decaying tail that is detected over 200 s 
from T . 

We simultaneously fit the spectral data of the Nal de- 
tectors with a source angle less than 60° (Nals 0, 1, 2, 3, 
and 5) and the data of the brightest BGO detector (BGO 
0) using the analysis package Rmfit 3.3rc8. An effective 
area correction is applied between each of the Nals and 
BGO during the fit process. This correction is used 
to handle possible discrepancies between the flux in the 
detectors due to the choice of the model to generate the 
instrument responses for instance. 
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Fig. 1. — The top two panels show the signal count rates as a 
function of time, as measured by the Fermi GBM detectors, from 
8 to 200 keV in Nal (top) and from 200 keV to 40 MeV in BGO 
(middle) . The bottom panel shows the evolution of the Band func- 
tion Epeajj (in blue) and the BB temperature kT (in red) over the 
duration of the burst. The vertical dashed lines indicate the period 
used in the time-integrated analysis. 

We performed a time-integrated spectral analysis over 
the main part of the burst (T -1.024s to T +83.969s) 
using the Band function. The Band parameters are 
in part fairly typical of the ensemble of GRBs, with 
a = -Q.67 ±0.01 and E pea k = 352 ± 6 keV (jPreece et all 
l2000t iKaneko et al.l I2006D . However, with an index 
/3 = —1.99 ± 0.02, the high-energy power law systemati- 
cally overshoots the observed flux above 1 MeV in BGO, 
as can be seen by the fit residuals in Figure [2] (top two 
panels), which also indicate systematic patterns at low 
energy. This suggests a simple Band function does not 
adequately represent the spectrum of this burst. 

We identify the best shape to fit the above-mentioned 
spectral deviation by fitting the same data simultane- 
ously with a Band function combined with each of the 
following models: single power- law (PL), Black Body 
(BB), Band function, power law with exponential cut off 
("Comp" for Comptonized model), and Gaussian. We 
select the best model by choosing the fit with the lowest 
Castor C-stat value (later C-stat). C-stat differs from 
the Poisson likelihood statistic by an offset which is a 
constant for a particular dataset. Table Q] shows the re- 
sults of these fits. The effective area correction described 
above is on the order of a few percent and does not change 
the C-stat for each fit more than a few units, nor does 
it change the value of the parameters resulting from the 
fit. 

While spectral deviations from the standard Band 
function were previously identified in the form of an ad- 
ditional PL to the Band function sometimes extending 
from the lo wer energy in the GBM to the higher ener gy 
in the LAT (lAbdo et al.ll2009t lAckermann et al.1l2010al lbi: 
iGuiriec et al.ll2010T ). in the case of GRB 100724B a PL 
spectral component does not improve on the Band-only 
fit and an additional BB component to the Band function 
is the best model to fit the spectral deviation. An equal 
C-stat is obtained for Band+Band and Band+Compt, 
but with a close to +1 for the additional Band and 
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Fig. 2. — The time-integrated spectrum of GRB 10072^B fit by a Band function (top two panels) and a Band + blackbody function 
(bottom bottom two panels). The left plots show the count spectra for the two models and the right plots the corresponding deconvolved 
vF v spectra. The data points appear as color crosses. Dashed lines indicate the individual spectral functions and solid lines shows the 
summed model fit. The addition of a BB spectral component over the brightest part of the burst (To-1.024s to To+83.969s) shows a 
significant improvement in the fit compared to a Band function by itself, particularly noticeable as the removal of trends with energy in 
the residuals compared to the Band-only fit. The region between 30 and 40 keV is excluded from the fit owing to calibration issues around 
the k-edge of the Nal detectors. We have verified this exclusion doesn't affect the recovered parameter values. 



TABLE 1 

Fit of the time-integrated spectrum of GRB 100724B from To + 1.024s to To+83.969s. The count spectrum using the NaI 

DETECTORS 0, 1, 2, 3, 5 AND BGO DETECTOR IS FIT SIMULTANEOUSLY WITH A STANDARD BAND FUNCTION AND WITH AN ADDITIONAL 
MODEL TO EVALUATE THE SHAPE OF THE SPECTRAL DEVIATION. BAND+BB IS PREFERRED OVER ALL THE OTHER COMBINATIONS. 



Models 


Standard Model 


Additional Model 




Band 


BB 


Compt 


Band 


Gaussian 


PL 


Cstat/dof 


Parameters 


Bpcak 


a 


P 


kT 


Bpcak 


index 


E 


a 





Centroid 


log 10 FWHM 


index 




Band 


352 
±6 


-0.67 
±0.01 


-1.99 
±0.01 




















1133/704 


Band+BB 


615 
±29 


-0.90 
±0.02 


-2.11 
±0.02 


38.14 
±0.87 


















1038/702 


Band+Compt 


708 
±48 


-0.94 
±0.02 


-2.13 
±0.02 




164 
±7 


+0.81 
±0.20 














1039/701 


Band+Band 


716 
±48 


-0.94 
±0.02 


-2.13 
±0.02 








60 
±7 


0.76 
±0.21 


< -5 








1039/700 


Band+Gaussian 


403 
±8 


-0.75 
±0.01 


-2.02 
±0.01 














103 
±2 


0.25 
±0.03 




1060/701 


Band+PL 


341 
±9 


-0.63 
±0.05 


-1.99 
±0.01 


















-1.93 
±1.59 


1131/702 
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Compt functions, and a very low value for j3 (only con- 
strained as an upper limit, below -5) for the extra Band 
function, the Band and Compt functions can be inter- 
preted as a Planck function. Even with more parame- 
ters, the additional Band and Compt functions resem- 
ble a BB component, reinforcing Band+BB as the best 
combination. Additio nal models were trie d, such as a 
log-Parabola function (lMassaro et al.ll2010() . but the re- 
sults were highly disfavored, and we exclude them from 
Table 1. 

Figure [2] (bottom two panels) shows the BB contribu- 
tion below Ep Ca k. Compared to the Band-only fit, E pea k 
is shifted towards higher energy to 615 ± 29 keV and 
f3 is lower with a value of —2.11 ± 0.02. This index is 
consistent with the flux de tected above 1 MeV , and the 
spectrum seen in the LAT (jTanaka et a,l.lf20Toh at higher 
energies, a is also significantly lowered to —0.90 ± 0.02. 

While the simultaneous fit of all the selected detectors 
provides the best constraints on the two spectral compo- 
nents, fits with Band+BB to combinations of individual 
Nal detectors with BGO result in similar parameter 
values and offer significant improvement over the Band- 
only fit. This provides a check that the BB component 
is real and not introduced by effects such as detector 
deadtime or spectral distortions that would affect each 
detector in a different way depending on the angle of the 
detector to the source. 

To verify that the improvement in the fit obtained by 
adding a BB component to the Band function is not 
a statistical fluctuation, we generated 20,000 synthetic 
spectra for each selected detector. For the simulations 
we used the parameters from the fit performed with the 
Band-only function, which we take as the null hypothe- 
sis. To create the simulated spectra, for each detector the 
real background is added to the source spectrum model 
and Poissonian fluctuations are applied to the sum. All 
the detectors are then fit simultaneously with both Band 
and Band+BB, and their C-stat are compared. None of 
the 20,000 simulated spectra give a difference larger than 
45 units of C-stat ([Band]-[Band+BB]), while in the real 
data, this difference is 95 units, corresponding to a prob- 
ability lower than 5 x 10~ 5 that the BB excess is due to 
statistical fluctuations. 

To study the evolution of the spectral components, 
22 time intervals were devised by requiring that each 
interval produce a Band+BB spectral fit with well- 
constrained Band function parameters, while attempt- 
ing to separate the peaks and valleys of the light curve 
so that the spectral fit parameters can be tracked with 
burst flux as well as with time. The bottom panel of 
Figure [1] exhibits the evolution of E poa k and BB tem- 
perature, kT, through these intervals. We notice that 
Epoak tracks the lightcurve and globally decreases over 
time. The BB component is detected throughout the 
burst, and its temperature shows weak correlation with 
Epeak- The significance of this correlation is difficult to 
assess, mostly because the variation in temperature is 
small. Overall, it appears the temperature is quite sta- 
ble, with Figure [3] showing more clearly the small scatter 
in kT. 

3. DISCUSSION 

With a fluence of ~ 5.2 x 10 -4 erg cm -2 measured 
in 85 sec from T -1.024s between 8 keV and 40 MeV, 
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Fig. 3. — Distribution of the time-resolved black-body tempera- 
ture kT. The spread in temperatures measured in the black-body 
component over the 22 time intervals shown in Figure [T] can be fit 
by a Gaussian distribution with a mean of 35.2 ± 2.3, and a lcr 
standard deviation of 6.0 ± 2.0. 

GRB 100724B is the most intense GRB detected by 
GBM over this energy range through 2010 September. 
Combined with the broad energy range of the GBM, this 
allows for accurate modeling of its energy spectrum even 
with this complex model. Previous observational results 
regarding thermal components in GRBs were ambigu- 
ous and some were limited to individual fine time-slices 
rather than a spectral fit over the entire emission pe- 
riod. Some studies showing BB fits did not demonstrate 
that the BB fit was stat istically preferred to a simple 
non-thermal component (|Ghirlanda et al.ll2003f ). Other 
analyses found BB+PL spectra for isolated portions of 
selected GRBs, raising the possibility that these spectra 
are actually adequately fit with a standard Band func- 
tion but that due to a weak signal in small time slices and 
extreme parameters for the Band func tion, a BB shape 
is competitive wi th the Band function ((Rvdc 2004, 120051 : 
IRvde et al.ll2r7I(l . The non-thermal component fit with 
a single power law suggested a break well beyond the 
common E pea k values, and the BB temperature and its 
variations intriguingly matched those of a typical E pca k. 
Despite the broader energy range of RHESSI GRB obser- 
vations, one analysis found difficult ies in fitting combined 
thermal plus non-thermal models (Bcllm 201Q). 

We find here that the joint BB plus non-thermal 
(Band) fit is highly statistically preferred so that in si- 
multaneously detecting both components we are con- 
fident in their correct identification. Time-resolved 
spectroscopy of GRB 100724B reveals that this BB 
component is seen throughout the burst and doesn't 
evolve much over time, while the n on-thermal com- 
ponent follows the typical variations ()Ford et al.l 119951 : 
iGuiriec et al.l I2010T ) . The consistency of the mean 
kT value with the temperature obtained in the time- 
integrated spectral fit, combined with the detection of 
the BB component throughout the burst, strengthen the 
case for an underlying thermal component in the gamma- 
ray emission seen from GRB 100724Q and show that the 
presence of the BB in the time-integrated spectrum can- 
not be attributed to spectral evolution of the Band func- 
tion during the burst. 
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Epcak varies substantially, from ~ 90 to ~ 1100 keV. At 
the same time, the thermal component remains relatively 
steady with the temperature varying only modestly be- 
tween 30 and 50 keV as suggested by a T oc L 1 / 2 depen- 
dence expected from a BB component. Time-averaged 
values of the temperature and the flux of the thermal 
component, and of the ratio of this flux over the total 
gamma-ray flux are 

kT = 38 ± 4 keV, F bb = (2.6 ± 1.4) x 10~ 7 erg/s/cm 2 
and F hh /F tot = 0.04 ± 0.02. 

In the standard fireball model these observables allow 
determination of the physical properties of the outflow 
and its photosphere. Owing to the imprecise and delayed 
localization of GRB 100724&, optical follow up to deter- 
mine the distance to the source was impossible. For this 
reason the temperature of the BB can be translated into 
a real source temperature only as a function of source 
distance. We assume in the following argument a typical 
redshift 2 = 1. We find that the Lorentz factor is T ~ 
325^ 1 / 4 / I ^ T ' , the photospheric radius is R p h ~ 5.6 x 

10 11 cm £ -3 / 4 /nx^ 4 an d the radius at the base of the flow 

is Rn ~ 1.2 x lO^cmf -1 /^ (jDaigne and Mochkovitchl 
l2002tlPe'er et al.ll2007| ). Here £ is a geometrical factor of 
order unity and /nt is the efficiency of the mechanism re- 
sponsible for the non-thermal emission. With an extreme 
efficiency /nt = 1, these estimates are in good agreement 
with the typical values expected in the fireball model. 
The dependence on redshift is not strong: at z=3 (resp. 
8), T ~ 645 (resp. 1290), R ph ~ 1.1 x 10 12 cm (resp. 
1.4xl0 12 cm), andi?o ^ l.lxl0 7 cm (resp. 6.9xl0 6 cm). 

Using more realistic values for the efficiency, the radius 
Rq is the most altered, with Rq ~ (3.6 — 40 km) for 
/nt — 0.1 — 0.5. Such small values are puzzling. If the 
central engine is a rotating black hole, as in the popu- 
lar collapsar model for long GRBs ()Wooslevlll993l) . with 
a minimal mass in the range 5-10 M@, such radii are 
smaller than the typical value expected for the innermost 
stable orbit, from 44-89 km for a non-rotating black hole 
to 22-43 km for a highly-rotating black hole having a 
spin a — 0.8. These results for the time-integrated spec- 
trum imply a small Rq or a very large efficiency and the 
constraint is even stronger in some time bins. We con- 
clude that observations of GRB 10072 ^B require either 
a very high efficiency for the non-thermal process, or a 
very small size of the region at the base of the flow, both 
of which are quite challenging for the standard fireball 
model, if not excluding it. 

A simple solution to this discrepancy between the stan- 
dard fireball model and the observations is to assume 
that the initial energy release by the central engine is 
not purely thermal, but that the flow is highly magne- 
tized close to the sourc e (jDaigne and Mochkovitchll2002l ; 
iZhang and FVer1[2009h . The magnetization a is the ratio 
of the Poynting flux over the power (thermal + kinetic) 
carried by the baryons. If no magnetic dissipation oc- 
curs below the photosphere, the efficiency /nt in the 
estimates of T, i? p h, and Rq above should be replaced by 
(1 + ct)/nt- A magnetization a > 1 will therefore rec- 
oncile the observations with physically acceptable values 
for the radius at the base of the flow and the efficiency of 
the mechanism responsible for the non-thermal emission. 
A similar conclusion is reached for scenarios where mag- 



netic dissipation occurs early and contributes efficiently 
to the acceleration of the jet. However, the appearance 
of a low intensity thermal component in the spectrum 
probably excludes the most extreme version of the mag- 
netized outflow scenario, where the energy is released by 
the central engine as a pure Poynting flux (a — oc). 

4. CONCLUSION 

We have shown that the simultaneous presence of 
thermal and non-thermal components to the spectra of 
GRB 100724B is statistically preferred. Although the 
non-thermal component is dominant, the black body flux 
is well within the GBM sensitivity. Deviations from the 
Band function may be measurable in less fluent bursts or 
in bursts where the thermal component is less prominent, 
providing that the black body component lies in the band 
pass of the instrument and its peak in energy is distin- 
guishable from -Epcak- If the presence of an unresolved 
thermal component in other bursts modifies the Band 
function parameters in the same sense as the Band-only 
fit for GRB 100724B, then we might expect a systematic 
bias yielding values of a and 8 that are higher (harder) 
than in the true non-thermal component. Two impor- 
tant consequences of this bias are that the perceived 
violation of the synchrotron limit that disallows values 
a > —2/3 (for slow-cooling electrons) and a > —3/2 (for 
fast-cooli ng electrons) may n ot be as common as sug - 
gested by iPreece et ah! (| 19981 ) and iCrider et ail ()1997fi , 
and that the relatively low rate of bur sts detected by the 
LAT compared to the predictions of IBand et al.l (2009) 
and the observations in Abdo et al. (in preparation) 
based on extrapolations of 8 from lower energies might 
be exp lained by this b i as in 8, a possibility suggested 
also bv lRvde and Pe'erl (f2009h . 

Our observations provide strong evidence for the pres- 
ence of a photospheric spectral component, long sus- 
pected to exist in the standard fireball model. In ad- 
dition, our results require implausible parameters for the 
standard baryonic fireball model and therefore favor a 
substantial magnetic component to the outflow. 

5. ACKNOWLEDGMENTS 

The GBM project is supported by the German Bun- 
desministerium fur Wirtschaft und Technologie (BMWi) 
via the Deutsches Zentrum fur Luft- und Raumfahrt 
(DLR) under the contract numbers 50 QV 0301 and 50 
OG 0502. 

A.J.v.d.H. was supported by an appointment to the 
NASA Postdoctoral Program at the MSFC, administered 
by Oak Ridge Associated Universities through a contract 
with NASA. 

S.F. acknowledges the support of the Irish Research 
Council for Science, Engineering and Technology, co- 
funded by Marie Curie Actions under FP7. 

P.M. acknowledges the support of NASA 
NNX08AL40G. 

F.R. acknowledges the support of the Swedish National 
Space Board. 



6 



Guiriec et al. 



REFERENCES 



Abdo ct al.: 2009, ApJ 706, L138 
Ackcrmann ct al.: 2010a, ArXiv c-prints 
Ackcrmann et al.: 2010b, "submitted to Apj" 
Band, D. et al.: 1993, ApJ 413, 281 
Band, D. et al.: 2009, ApJ 701, 1673 
Bcllm, E. C: 2010, ApJ 714, 881 

Bhat, P. N.: 2010, GRB Coordinates Network, Circular Service, 

10977, 1 (2010) 977, 1 
Crider, A. ct al.: 1997, ApJ 479, L39+ 

Daigne, F. and Mochkovitch, R.: 2002, MNRAS 336, 1271 
Ford, L. A. et al.: 1995, ApJ 439, 307 
Gehrels, N.: 1997, Nuovo Cimcnto B Scric 112, 11 
Ghirlanda, C, Cclotti, A., and Ghisellini, C: 2003, AfcA 406, 
879 

Golcnctskii, S. ct al.: 2010, GRB Coordinates Network, Circular 

Service, 10981, 1 (2010) 98171 
Goodman, J.: 1986, ApJ 308, L47 
Guiriec, S. et al.: 2010, ApJ 725, L73-L78 
Kancko, Y., Preece, R. D., Briggs, M. S., Paciesas, W. S., 

Meegan, C. A., and Band, D. L.: 2006, ApJS 166, 298 
Massaro, F., Grindlay, J. E., and Paggi, A.: 2010, ApJ 714, L299 
Meegan, C. ct al.: 2009, ApJ 702 , 791 
Meszaros, P.: 2002, ARAfcA 40, 137 
Meszaros, P. and Rccs, M. J.: 2000, ApJ 530, 292 



Pe'er, A., Ryde, F., Wijers, R. A. M. J., Meszaros, P., and Rees, 

M. J.: 2007, ApJ 664, LI 
Piran, T.: 1999, Phys. Rep. 314, 575 
Piran, T.: 2004, Reviews of Modern Physics 76, 1143 
Preece, R. D., Briggs, M. S., Mallozzi, R. S., Pendleton, G. N., 

Paciesas, W. S., and Band, D. L.: 1998, ApJ 506, L23 
Preece, R. D., Briggs, M. S., Mallozzi, R. S., Pendleton, G. N., 

Paciesas, W. S., and Band, D. L.: 2000, ApJS 126, 19 
Rees, M. J. and Meszaros, P.: 2005, ApJ 628, 847 
Ryde, F.: 2004, ApJ 614, 827 
Ryde, F.: 2005, XpJ 625, L95 

Ryde, F., Axelsson, M., Zhang, B. B., McGlynn, S., Pe'er, A., 
Lundman, C, Larsson, S., Battelino, M., Zhang, B., Bissaldi, 
E., Bregeon, J., Briggs, M. S., Chiang, J., de Palma, F., 
Guiriec, S., Larsson, J., Longo, F., McBreen, S., Omodei, N., 
Petrosian, V., Preece, R., and van der Horst, A. J.: 2010, ApJ 
709, L172 

Ryde, F. and Pe'er, A.: 2009, ApJ 702, 1211 

Tanaka, Y., Ohno, M., Takahashi, H., Uchara, T., Omodei, N., 

Chiang, J., and Guiriec, S.: 2010, GRB Coordinates Network, 

Circular Service, 10978, 1 (2010) 97871 
Woosley, S. E.: 1993, ApJ 405, 273 
Zhang, B. and Pe'er, A.: 2009, ApJ 700, L65 



